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The pion exchange between nucleons generates a strong tensor interaction, which pro- 
vides a large attractive contribution for the binding energy of nucleus. This non central 
tensor interaction is difficult to handle in the shell model framework, which hinders full un- 
derstanding of nuclear structure. We develop the tensor-optimized shell model (TOSM) for 
the strong tensor interaction and now we are able to use bare nucleon-nucleon interaction 
with the help of the unitary correlation operator method (UCOM) for the short-range hard 
core. We adopt the nucleon-nucleon interaction, AV8', and calculate explicitly the ground 
state of 4 He and make a detailed comparison with rigorous few-body model calculations. We 
show a large amount of success of the tensor-optimized shell model with bare nucleon-nucleon 
interaction for 4 He. 

(N 

§1. Introduction 

00 

. It is important to develop a theoretical framework to calculate nuclear struc- 

ture with many nucleons using the realistic nucleon-nucleon interaction, which is 
obtained from two nucleon scattering. Recently, it has become possible to calculate 
nuclei up to a mass of approximately A ~ 12 1 ^ 3 ^ using the realistic nucleon-nucleon 
interaction. The method used for the calculation is the Green's function Monte- 



Carlo method (GFMC) with the use of relative nucleon coordinates. This method 
introduces various correlation functions with many variational parameters in the nu- 
clear wave function. In GFMC, the nuclear structures and binding energies were 
successfully reproduced by including three-body interaction. One big surprise is the 
extremely large contribution of the one pion exchange interaction, which is about 70 
~ 80% of the entire nucleon-nucleon interaction. In principle, they can extend this 
method to calculate heavier nuclei. It is, however, extremely time-consuming even 
with the present computer power. Hence, it is strongly desired to develop a new 
method of calculating nuclei with large nucleon numbers using the nucleon-nucleon 
interaction. 

The nucleon-nucleon interaction has distinctive features, namely there exist 
strong tensor interaction at intermediate distance caused by pion exchange and 
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strong short-range repulsive interaction at a short distance caused by quark dy- 
namics. Although these two interactions have totally different characteristics, it is 
customary to adopt the Brueckner Hartree-Fock theory to integrate out the high- 
momentum components on the same footing and use the resulting G-matrix as an 
effective interaction in the shell model. In this way, we lose information on the tensor 
correlation and a short-range correlation in the shell model wave function. Hence, we 
search for a powerful method of treating treat explicitly both the tensor interaction 
and the short-range interaction to study of not only light nuclei but also medium 
and heavy nuclei. 

There have been two important developments for this purpose. One is the finding 
that the tensor interaction is of intermediate range, and hence, we can express the 
tensor correlation in a reasonable shell model space. 4 )' 5 ) We call this method Tensor- 
Optimized Shell Model (TOSM), wherein the nuclear wave function is written in 
terms of the standard shell model state and a sufficient amount of two-particle two- 
hole (2p2h) states. This TOSM formalism is based on the success of the parity and 
charge projection in the treatment of the pion exchange interaction. 6 )' 7 ) We have 
shown that the tensor interaction could be treated properly by taking a reasonable 
amount of multipoles (/ < 5) in the 2p2h wave functions with the optimization of 
the radial parts of the particle states. The other is the Unitary Correlation Operator 
Method (UCOM) for the treatment of the short-range correlation. 8 )' 9 ) The short- 
range repulsive interaction is of very short-range and it is suited to treat the short- 
range correlation using unitary transformation and take the approximation to use 
only up to the two-body operators. This approximation is justified because the 
volume associated with the short-range correlation is extremely small, where more 
than three nucleons rarely enter the small volume. This is not the case for the tensor 
correlation, since the tensor interaction is of intermediate and long range as discussed 
by Neff and Feldmeier. 9 ) 

Our idea is to combine these two methods, TOSM and UCOM, to develop a 
theoretical framework that can describe medium and heavy nuclei beyond the light 
nuclei using the realistic nucleon-nucleon interaction. We can use the TOSM for the 
strong tensor interaction utilizing the intermediate nature caused by finite angular 
momentum of the relative wave function and the UCOM for the strong short-range 
interaction utilizing the short-range nature. We use completely different methods 
for these two distinctive characteristics of the nucleon-nucleon interaction. After 
demonstrating its power, we hope to apply the newly developed method, which we 
call TOSCOM, to many nuclei. Using TOSCOM, we aim to understand the roles of 
the tensor and short-range correlations in nuclei using bare interaction. As a good 
start, we would like to apply TOSCOM to 4 He. Hence, there are two purposes of this 
study. One is to see how this method works for the treatment of the bare nucleon- 
nucleon interaction. The other is to compare the obtained results with rigorous 
calculations to check the accuracy of TOSCOM. From this comparison, we can see 
how far we can describe the short-range and tensor correlations and to determine 
what we need to do for the further improvement of TOSCOM in order to solve the 
nucleus as precisely as possible. 

There are several methods for the description of few-body systems. These 
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methods are compared each other for 4 He using the AV8' nucleon-nucleon inter- 
action. 3 )> 10 ) These methods are called the Faddeev-Yakubovsky equations (FY), 11 ) 
the coupled-rearrangement-channel Gaussian-basis variational method (CRCVM), 12 ) 
the stochastic variational method(SVM), 13 ) the hyperspherical harmonic method 
(HH), 14 ) GFMC, 15 ) the no-core shell model (NCSM), 16 ) the effective interaction hy- 
perspherical harmonic method (EIHH). 17 ) Although the methods of numerical calcu- 
lations are largely different, these methods provide essentially the same results for the 
4 He structure and also the amount of kinetic energy and components of the nucleon- 
nucleon interaction. All the methods described in Ref. 10) have a common feature 
where the wave functions are expressed in terms of the relative nucleon coordinates 
(Jacobi coordinate). We shall call this method as described here the T-coordinate 
method, since the relative coordinate between two nucleons is the coordinate of the 
nucleon-nucleon interaction. Hence, it is easy to describe the correlations between 
two nucleons by taking sufficient variational variables. This T-coordinate method is 
advantageous for rigorous calculations when there exists sufficient computer power. 

On the other hand, we would like to develop a theoretical framework to describe 
wave functions in terms of single-particle coordinates, which we call U-coordinate 
method. This ^-coordinate method can be used to describe nuclei with many nu- 
cleons relatively easier than the T-coordinate method. Furthermore, we are able 
to describe the wave function on the basis of the shell model picture, and hence, 
it becomes easier to interpret the calculated results in the shell model sense. The 
difficulty, on the other hand, is to express the correlations of the relative motion 
between two nucleons, which are caused by the short-range repulsive interaction and 
the tensor interaction in nucleon-nucleon interaction. We overcome this problem by 
developing TOSCOM to describe the short-range and tensor correlations simultane- 
ously. 

We also mention the differences of TOSCOM from the other works based on 
the shell model. Recently, NCSM result has been reported for 40 Ca starting from 
the realistic interaction, and discussions are made for the convergence of the binding 
energy 18 ^ In their scheme, transformed interactions such as V\ ow k, using UCOM 
and Lee-Suzuki technique, 16 ) are employed. In these interactions, high-momentum 
components of the short-range and tensor correlations are renormalized. Thus, the 
truncated shell model space calculation is applicable to discuss the binding energy. 
On the other hand, in TOSCOM, the truncation of the model space is not introduced. 
For the tensor interaction, we can directly use the bare tensor interaction to evaluate 
the matrix elements and discuss the characteristics of the tensor interaction in the 
nuclear structure explicitly. For the short-range part, we employ UCOM, in which 
we can also obtain the explicit wave function including the short-range correlation 
through the UCOM transformation. 

Otsuka et al. investigated the role of the 'tensor interaction', in particular, for 
neutron-rich nuclei. 19 )> 20 ) They calculated the tensor interaction matrix elements 
mainly arising from the exchange term of the tensor interaction in the Hartree-Fock 
scheme. On the other hand, we go beyond the Hartree-Fock approximation and treat 
the dominant part of the tensor interaction by taking the 2p2h wave functions up to 
all orders for the particle states. In particular, the 0p0h-2p2h coupling of the tensor 



4 



T. Myo, H. Toki and K. Ikeda 



interaction is essential for describing the dominant part of the tensor interaction, 
which provides about half of the interaction matrix element. In TOSCOM, this 
2p2h contribution is explicitly included, and we focus on the strong tensor correlation 
represented by the 2p2h wave functions. 

This paper is arranged as follows. In 52 and we describe the TOSM + 
UCOM (TOSCOM) in detail. We shall apply TOSCOM to 4 He in g] and compare 
the numerical results with those of the few-body methods. In this study is 
summarized and perspectives are described. 

§2. Tensor-optimized shell model 

We shall begin with many-body Hamiltonian, 

B = J2 T i- T cm+J2 V *3 ^ 

i i<j 

with 

Vij=v C + v T +v LS + v Clm b _ (2 . 2) 

Here, Tj is the kinetic energy of all the nucleons with T cm being the center of mass 
kinetic energy. We take the bare nucleon-nucleon interaction for Vij such as the AV8' 
consisting of central (vfj), tensor (vj^) and spin-orbit {v^ s ) terms. The v^ mh is the 
Coulomb term. We describe the many-body system with many-body wave function, 
by solving the equation H\P = E'F. In TOSCOM, we take the ^-coordinates 
to express \P '. We explain here the typical features of nucleon-nucleon interaction, 
tensor interaction and short-range repulsion, and how we treat them in TOSCOM. 
The pion exchange interaction is a long- and intermediate-range interaction, which 
contains strong tensor interaction. This tensor interaction has a large strength in 
the intermediate range. 4 ) Hence, we hope to describe the tensor correlation in terms 
of a reasonable amount of multipoles of single-particle states in the 2p2h states with 
high-momentum component. In fact, we have shown that the tensor correlation is 
expressed in terms of multipoles up to I < 5 for 4 He. 4 ) On the other hand, there 
exists a strong repulsive interaction in the short-range part of the nucleon-nucleon 
interaction. This is the other difficulty to be considered for the nuclear many-body 
problem. For this problem, Feldmeier et al. have demonstrated that UCOM can be 
used to treat the short-range correlation. 8 )' 9 ) 

We begin with TOSM and write the case of 4 He explicitly as an example. The 
wave function, is written as 

V = C \0)+Y / C P \2 P 2h) p . (2-3) 

p 

Here, the wave function |0) is a shell model wave function and | (Os) 4 ) for 4 He. \2p2h) 
represents a 2p2h state with various ranges for the radial wave functions of particle 
states. We can write \2p2h) as 

\2 P 2h) p = \mi(xim(x2)] J [C^i)C 4 4 (x 2 )] J ]°>A • (2-4) 
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The suffix A of the wave function indicates anti-symmetrization of the wave functions. 
Here, p denotes representable quantum number of 2p2h states, which are expressed 
with particle (hole) wave functions ip™ {ip^)- The index, n, is to distinguish the 
different radial components of the single-particle wave function, tp. The index, a, 
is a set of three quantum numbers, j and t z , to distinguish the single-particle 
orbits, where I and j are the orbital and total angular momenta of the single-particle 
states, respectively, and t z is the projection of the nucleon isospin. The normalization 
factors of the two particle states are included in the wave functions given in Eq. (|2-4p . 
For 4 He, the coupled spin, J, of two nucleons is J = or J = 1. We omit writing 
the coupled isospin, which should be either or 1 depending on the value of J. We 
have used Gaussian functions for radial wave functions to express more effectively 
compressed radial wave functions. 4 ) The shell model technique is used to calculate 
all the necessary matrix elements, which are expressed explicitly in the Appendix. 

We explain the Gaussian expansion technique for single-particle orbits. 12 )' 21 ) 
Each Gaussian basis function has the form of a nodeless harmonic oscillator wave 
function (HOWF), except for Is orbit. When we superpose a sufficient number 
of Gaussian bases with appropriate length parameters, we can fully optimize the 
radial component of every orbit of every configuration with respect to the total 
Hamiltonian in Eq. f)2 • 1 j) . We construct the following ortho-normalized single-particle 
wave function ip™ with a linear combination of Gaussian bases {4> a } with length 
parameter b a ^ m . 

N a 

V<» = Yl <m Mr, b a , m ) for n = 1, ■ ■ ■ , N a (2-5) 

m=l 

Here, N a is the number of basis functions for a, and m is an index that distinguishes 
the bases with different values of b a ^ m . The explicit form of the Gaussian basis 
function is expressed as 

Mr,b a , m ) = iV,(6«,m) r l e~^ b ^ 2 / 2 \Yi(r) t ^ /2 ] jX t., (2-6) 



Ni{b a , 



o z,-(2*+3) 
4 Ua,m 

r{i + 3/2) 



(2-7) 



The coefficients {d™ m } are determined by solving the eigenvalue problem for the 
norm matrix of the non orthogonal Gaussian basis set in Eq. (|2-6p with the dimension 
N a . Following this procedure, we obtain new single-particle wave functions {V'q} 
using Eq. (|2-5p . 

We choose the Gaussian bases for the particle states to be orthogonal to the 
occupied single-particle states, which is 0si/ 2 in the 4 He case. For 0si/ 2 states, we 
employ one Gaussian basis function, namely, HOWF with length 6os 1/2 ,m=i = bo s . 
For lsi/2 states, we introduce an extended Is basis function orthogonal to the 0si/ 2 
states and possessing a length parameter 6i Sjm that differs from 6os- 4 ^ In the extended 
Is basis functions, we change the polynomial part from the usual Is basis states to 
satisfy the conditions of the normalization and the orthogonality to the Os state. 
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Two-body matrix elements in the Hamiltonian are analytically calculated using 
the Gaussian bases, 12 )' 21 ) whose explicit forms are given in the Appendix for central, 
LS and tensor interactions, respectively. In the numerical calculation, we prepare 9 
Gaussian functions at most with parameters of various ranges to obtain a convergence 
of the energy 

Furthermore, we have to take care of the center-of-mass excitations. For this 
purpose, we use the well-tested method of introducing a center-of-mass term in the 
many-body Hamiltonian. 22 )> 23 ) 



H — x (^ + 5 WJ ^ -§*")■ (2 - 8) 

A 1 A ^ 

Pern = J>>, ^ = ^ W = ( 2 - 9 ) 
i=l i=l ® s 



Here, m and A are the nucleon mass and the mass number, respectively, and bo s is 
the length parameter of HOWF for the hole Os state. We take a sufficiently large 
coefficient, A, to project out only the lowest HO state for the center-of-mass motion. 
In the numerical calculation, the excitation of the spurious center-of-mass motion is 
suppressed to be less than 10 keV. 

The variation of the energy expectation value with respect to the total wave 
function •P'( 4 He) is given by 



= ' (2 ' 10) 



which leads to the following equations: 

d(&\H-E\&) n d(¥\H-E\V) 



dba^m dC p 



(2-11) 



Here, E is a Lagrange multiplier corresponding to the total energy. The parameters 
{b a ,m} for the Gaussian bases appear in non linear forms in the energy expectation 
value. We solve two types of variational equations in the following steps. First, 
fixing all the length parameters b a>m , we solve the linear equation for {C p } as an 
eigenvalue problem for H with partial waves up to L max . We thereby obtain the 
eigenvalue E, which is a function of {6 a>m }. Next, we try to search various sets of 
the length parameters {6 a ,m.} to find the solution that minimizes the total energy. In 
this wave function, we can describe the spatial shrinkage with an appropriate radial 
form, which is important for the tensor correlation. 4 ) 
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§3. Formulation of UCOM 

3.1. short-range correlation in UCOM 

We employ UCOM for the short-range correlation. Feldmeier et al. worked out 
a unitary correlation operator in the form, 8 )' 9 ) 

C = exp I -i 9ij J = Yl aj (3-1) 

with Cij = exp(— i gij). Here, i and j are the indices to distinguish particles. Here, 
the two-body operator, g^, is a Hermite operator, and hence, C is a unitary operator. 
We express the full wave function, in terms of less sophisticated wave function, <P, 
as ^ = C<&. Hence, the Schroedinger equation, = becomes H<P = E<&, where 
H = C^HC. If we choose properly the unitary correlator, C, we are able to solve 
more easily the Schroedinger equation. Moreover, once we obtain <P, we can then 
obtain the full wave function, by the unitary transformation, ^ = C<P. Since C is 
expressed with a two-body operator in the exponential, it is a many-body operator. 
In the case of the short-range correlation, we are able to truncate modified operators 
at the level of two-body operators. 8 ) 

In the actual calculation of UCOM, we define the operator as 

9ij = \ {Pr,ijs(rij) + sir^pr^} , (3-2) 

where the momentum p r ^j is the radial component of the relative momentum, which 
is conjugate to the relative coordinate r^-. s{rij) is the amount of the shift of the 
relative wave function at the relative coordinate, r^-, for each nucleon pair. Hereafter, 
we omit the indices i and j for simplicity. We also introduce R+(r) as 



s(0 ' 

which leads to the following relation, 

dR + (r) _ s(R + (r)) 
dr s(r) 



R+(r) Jt 

-1, (3-3) 



(3-4) 



In UCOM, we use R+{r) instead of s{r) to use the UCOM prescription. R+{r) 
represents the correlation function to reduce the amplitude of the short-range part 
of the relative wave function in nuclei and can be determined for four spin-isospin 
channels independently. The explicit form of the transformation of the operator for 
the relative motion is given as 

c t rc = R + (r), c f p r c = — | p r — | , cHc = I, (3-5) 
c f sc = s, c^ S12C = S12, c^v(r)c = v(R+(r)), (3-6) 
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where the operators I, s and S12 are the relative orbital angular momentum operator, 
intrinsic spin operator and the tensor operator, respectively. v(r) is the arbitrary 
function depending on r, such as potential. 

In the calculation using UCOM, we parametrize -R+(r) in the same manner as 
proposed by Neff-Feldmeier and Roth et al. 8 )' 9 )' 24 ) We assume the following forms 
for even and odd channels, respectively. 

R»r(r) =r + a(^y exp[- exp(r//3)], (3-7) 

#° dd (r) = r + a (1 - exp(-r/ 7 )) exp[- exp(r//?)] (3-8) 

Here, a, (3, 7 are the variational parameters to optimize the function R+(r) and 
minimize the energy of the system. They are independently determined for four 
channels of the spin-isospin pair. In the actual procedure of the variation, once 
we fix the parameters included in R + (r), we solve the eigenvalue problem of the 
Hamiltonian using Eq. (|2-lip and determine the configuration mixing of the shell 
model-type bases. Next, we try to search various sets of the R+(r) parameters to 
minimize the obtained energy. 

3.2. Extension of UCOM S-wave UCOM- 

In the framework of UCOM, we introduce the UCOM function R+(r) for each 
spin-isospin channel and ignore the partial wave dependence of R+(r). It is generally 
possible to introduce the partial wave dependence in UCOM and then R+(r) func- 
tions are determined in each relative partial wave in the two-body matrix elements. 
Here, we consider the specific case of this extension of UCOM by taking care of 
the characteristics of the short-range correlation. One of the simplest cases of this 
extension is UCOM for only s-wave relative motion, since all the other partial waves 
I except for s-wave (I = 0) have r l behavior near the origin, where the short-range 
hard core is extremely large. Hence, this r l behavior largely cuts down the effect of 
the short-range hard core. However, only the s-wave function is finite at the origin, 
and the behavior in the origin is determined by the hard core dynamics. In fact, the 
method used by Feldmeier et al. is to determine the unitary operator to reproduce 
the short-range behavior of s-wave relative wave function. 

When we incorporate S-wave UCOM (S-UCOM, hereafter) into TOSM, we ex- 
tract the relative s-wave component in all the two-body matrix elements in TOSM 
using the U-type basis expanded by the Gaussian functions. For numerical calcula- 
tions, we prepare the completeness relation consisting of the T-type basis functions 
|T) as 

l = Y j m)(T i \, |T> = \WiJ£\v,xs]jXt), (3-9) 

i 

where the T-type basis is expanded by the two coordinates of the relative part r and 
the center of mass part R of two nucleons, which are the set of Jacobi coordinate. 
The orbital angular momenta of each coordinate, r and R, are I and L, respectively. 
It is easy to prepare the s-wave relative part by considering I as zero in the T-type 
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basis. We construct the above completeness relation of the T-type basis states by 
diagonalizing the norm matrix expanded by the finite number of Gaussian basis 
functions for two coordinates. In the actual calculation, we use 12 bases for each 
coordinate, with which convergence is achieved. 

We calculate the matrix elements of the arbitrary two-body operator O including 
the S-UCOM correlator C s using the V-type basis with indices a and (3. Here, we 
insert the above T type completeness relation in Eq. (|3-9p as 

(V a \C$OC s \V p ) = ^(V a \T^ ■ miCtOCslTj) ■ (TjlVp). (3-10) 

ij 

The matrix element using T-type base, {Ti\C\OC ' s \Tj) , is calculated for the two-body 
kinetic part and central, tensor interactions. For the kinetic part and the central 
interaction, the matrix elements conserve the relative angular momentum, and then 
we can easily calculate the matrix elements of the transformed operator C\OC s . 
For the tensor interaction, the sd coupling matrix elements are properly treated, in 
which C s is operated on only the s-wave part of the relative motion. In this case, 
the operator C s acts on the s-wave relative Gaussian basis function (j>i=o(r), which 
is transformed as 

C S &=o(r) = J=o (i2_(r)), (3-11) 

where R-(r) is the inverse transformation of R+(r), namely, R-(R+(r)) = r. The 
matrix elements of the T-type basis function are calculated using the above trans- 
formed wave function. We also calculate the overlap between V type and T-type 
bases using Gaussian basis functions, whose explicit form is given in the Appendix. 

§4. Numerical results of TOSCOM 

In this section, we show our numerical results for 4 He using AV8' potential, 
which consists of central, LS and tensor terms and is used in the rigorous calculation 
given by Kamada et al. where the Coulomb term is ignored. 10 ) 

4.1. Optimization of R+(r) 

First, we determine the UCOM functions R+(r) for the calculation of TOSCOM. 
In UCOM, we optimize R+(r) function by changing the three parameters of a, (3 
and 7 to search for the energy minimum in TOSCOM. In Table [U the optimized 
three parameters in S'-UCOM are listed. The demonstration of the calculated result 
to search for the energy minimum is shown in Fig. [TJ In Fig. [TJ the total energy 
is plotted as a function of the range parameter, (3 of R+(r) function of the triplet 
even channel in the case of L max being 10. We have already minimized a and 7 for 
each (3 in this calculation. We see clearly a desired behavior, where the energy has 
a minimum as a function of (3. Hence, we fix (3 that provides the lowest energy. In 
Fig. [21 R+{r) functions used in the present study are plotted in comparison with 
the case in Ref. 24). For the odd channel, in accordance with the discussion in Refs. 
9) and 24), we cannot find the optimum value of R+(r) in the two-body cluster 
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Table I. Optimized parameters in R+(r) in TOSCOM for four channels in fm in the present work. 





a 


P 


7 


singlet even 


1.32 


0.88 


0.36 


triplet even 


1.33 


0.93 


0.41 


singlet odd 


1.57 


1.26 


0.73 


triplet odd 


1.18 


1.39 


0.53 



approximation of the UCOM transformation for the Hamiltonian. Hence, we decide 
to fix the range of R + (r), namely, (3 as the same one adopted in Ref. 24) and optimize 
a and 7, while the variation of R+(r) for the odd channel does not have significant 
effects on the energy and other properties of 4 He in comparison with the original 
case. 9 )' 24 ) Essentially, two types of parameter set of R+(r) in the present study and 
Ref. 24) give the similar form of R+(r) for even channels, in which we omitted the 
correlation function for the even channels except for s- waves. This result indicates 
that the correlation functions for the short-range repulsion are uniquely determined 
for each channel. 

-20.0 
^ -20.5 

I - 21 -° 
>, -21.5 

o -22.0 

W -22.5 

-23.0 

0.8 1.0 1.2 1.4 1.6 

Pte [fm] 

Fig. 1. Energy surface of 4 He with respect to /3 of triplet even channel with L m ax being 10. 





r [fm] r [fm] 

Fig. 2. Short-range correlation functions, R + (r), for UCOM in even and odd channels. The thin 
curves are for the R+ (r) function of Roth et al. and the thick curves are for that of this work. 
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Fig. 3. Energy of 4 He in TOSCOM as a function of the maximum angular momentum, L max . 
The circles are the results where UCOM is used for all the partial waves. The triangles (S- 
UCOM(ap)) are the results using S-UCOM with the approximation neglecting the UCOM 
effect in the tensor interaction matrix elements. The squares are the results using S-UCOM 
accurately. 



4.2. TOSCOM results for 4 He 

Next, we show the calculated results of the energy of 4 He as a function of L max in 
Fig. [3l We shall then compare the obtained results with the benchmark calculation 
given in Ref. 10). To start with, we show the ordinary UCOM case where UCOM 
is used for all the partial waves. The calculated results of the energy are indicated 
in Fig. [3] by circles as a function of the maximum angular momentum, L max - The 
results show good convergence to reach —19 MeV, while the exact value of the few- 
body calculations is approximately —26 MeV as indicated in Fig. [3j We would like 
to point out that we can calculate the binding energy directly using the nucleon- 
nucleon interaction in TOSCOM. However, the binding energy is small. The tensor 
interaction matrix element is approximately —50 MeV. On the other hand, in the 
previous study, 4 ) we obtained approximately —60 MeV for the tensor interaction 
matrix element to check the validity of TOSM, when we used G-matrix for the 
central interaction to renormalize the short-range repulsion and retained the bare 
tensor interaction of AV8' in our previous calculation. This fact indicates that the 
treatment of the short-range repulsive interaction is interfering with the contribution 
of the tensor interaction. This is due to a large removal of the short-range part of the 
relative wave functions in UCOM, in particular, in the d-wave part of the sd coupling 
of the tensor interaction matrix element, where the tensor interaction possesses some 
amount of strength. We have also calculated the contributions beyond the 2p2h 
configurations in TOSM such as 3p3h and ApAh configurations. When we include 
the ApAh configurations within the p-shell, their contribution to the binding energy is 
approximately 50 keV. This fact denotes that these more complicated wave functions 
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contribute very little in the total 4 He wave function. 

We have decided to restrict the use of UCOM to the relative s-wave only for 
the even channel (S-UCOM), where the treatment of the short-range repulsion is 
absolutely necessary. In other partial waves, we have the centrifugal potential that 
cuts out the short-range part from the wave functions of the higher partial waves. 
In this case, we can use the modified interaction and the kinetic energy only for the 
relative s-wave component in the even channels. Since the use of the UCOM for the 
odd partial wave is slightly better, we use the UCOM for all odd partial waves. As a 
starting calculation, we have neglected the S-UCOM correlation in the calculation of 
the tensor interaction matrix elements. The results in this case (S'-UCOM(ap)) are 
indicated also in Fig.[3]by triangles. The energy converges to —24 MeV, which is now 
very close to the exact one as shown in Fig. [3l In this case, the tensor interaction 
matrix element is —61 MeV, which becomes close to the exact value of —68 MeV. 
This improvement mainly comes from the increase in the sd coupling of the tensor 
interaction matrix element, however, this calculation is still not yet perfect. We 
have to treat the effects of the short-range repulsion on the tensor interaction matrix 
element. Hence, we have worked out the formulation to treat the rigorous s-wave 
function with the effect of the short-range repulsion for the calculation of the tensor 
interaction matrix elements as explained in the previous section. 

The numerical calculation is quite involved in the 5-UCOM case, since the s- 
wave relative wave function with the effect of short-range repulsion should be used 
for the tensor interaction matrix element. We show the calculated results for the 
total energy by the squares in Fig. We see quite a satisfactory result for the 
total energy, which is approximately —22 MeV. We show now all the components 
of the energy for 4 He in Fig. HI All the energy components show the saturation 
behavior as function of L max . In the tensor component, the saturation is obtained 
at around L max being 8. For the other components, their saturation points are seen 
at the similar L max . A very interesting feature is the kinetic energy, which goes up 
to a large value as the tensor interaction matrix element becomes large. As for the 
comparison with the rigorous calculation, we see that V c satisfies the rigorous value, 
which is approximately —55 MeV. On the other hand, the tensor interaction matrix 
element, Vt converges to —55 MeV, while the rigorous one is —68 MeV. The kinetic 
energy is approximately 90 MeV, while the rigorous one is 102 MeV. The LS matrix 
element is also smaller than the rigorous value. As the net value, the total energy, E, 
is —22 MeV and the rigorous value is —26 MeV. A detailed comparison is performed 
in Table im in which the converged energies in TOSCOM are shown with the rigorous 
calculations. One of the possibilities for the lack of the energy in TOSCOM is due 
to the separate treatment of the short-range and tensor correlations. Although the 
dominant part of the tensor interaction is of intermediate and long range, there 
may remain some small strength in the short-range part of the tensor interaction, 
which can couple with the short-range correlations. This effect can be included by 
extending the truncation of the UCOM transformation in the Hamiltonian to more 
than the two-body level. Three-body term of the UCOM transformation is one of 
the possibilities to overcome the lack of energy in TOSCOM. 8 ) 

We show the properties of the 4 He wave functions obtained in the present cal- 
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Fig. 4. Matrix elements of the central interaction (Vc), tensor interaction (Vr) and the spin-orbit 
interaction (Vls) together with the kinetic energy (Kinetic) and total energy (Energy) in the 
Hamiltonian for 4 He as function of L max . We observe good convergence for all the matrix 
elements. These values are compared with the benchmark results of Ref. 10), which are indicated 
by the thick short solid lines on the right-hand side of the figure. 

Table II. Total energy, matrix elements of the Hamiltonian and radius of 4 He compared with the 
benchmark results in Ref. 10). Units are in MeV for the total energy and the matrix elements, 
and fm for the radius of 4 He. 





Energy 


Kinetic 


Central 


Tensor 


LS 


Radius 


Present (UCOM) 


-19.46 


88.64 


-56.81 


-50.05 


-1.24 


1.555 


Present(5-UCOM) 


-22.30 


90.50 


-55.71 


-54.55 


-2.53 


1.546 


FY in Ref. 10) 


-25.94 


102.39 


-55.26 


-68.35 


-4.72 


1.485 



Table III. Probabilities of the total intrinsic spin 5 components in %. 

= S = 5=1 5 = 2 
89.41 2.47 8.12 



culation in Table IIIII The three probabilities of the total intrinsic spin S are shown 
in this table, in which the coupled value of the intrinsic spin of four nucleons is cal- 
culated without the UCOM transformation. It is found that the S = 2 component 
is larger than the S = 1 case. This S = 2 component is strongly caused by the 
tensor interaction, which can change the spin of two nucleon pairs by two. We show 
also the expectation values of the potentials for four channels in Table [TVl The large 
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Table IV. Expectation values of the potentials for four channels. 





1 E Z E 1 


3 o 


Central 


-37.49 -19.13 0.47 


0.45 


Tensor 


-54.03 


-0.52 


LS 


-3.07 


0.54 



Table V. Mixing probabilities in the 4 He ground state in %. 





(Osfoo 


82.48 




(0s)- 2 (0p 1/2 ) 2 10 


2.54 


(0s 


r 2 [(ls 1/2 )(0d 3 /2)]io 


2.34 


(0s 


W [(0P3/2)(0/ 5 / 2 )]iO 
w[(0Pl/ 2 )(0p3/2)]l0 
ro 2 [(0rf5/2)(0 57 /2)]l0 
(0 S )r o 2 (0d 3 /2)lO 


1.90 


(0s 


1.55 


(0s 


0.79 




0.44 




remaining part 


7.96 



contributions of the triplet even channel ( 3 E) are shown for the tensor and LS terms. 
On the other hand, the odd channel contributions are very small. 

In Table El we list the mixing probabilities of the dominant configurations in 
4 He. The subscripts 00 and 10 represent J and T, the spin and isospin quantum 
numbers, respectively. It is found that the 2p2h configurations with (J, T)=(l, 0) for 
the particle pair state are significantly mixed. These spin and isospin are the same 
as those for the deuteron, and thus, this two-nucleon coupling can be understood as 
a deuteron-like correlation. 4 ) 

In Table ED we list the dominant components of the matrix element of the tensor 
interaction between the OpOh and 2p2h states. We can expand the total contribution 
of the tensor interaction {Vr) into two types of the matrix element between OpOh 
and 2p2h states and between 2p2h and 2p2h states of the wave function. It is found 
that the former 0p0h-2p2h coupling produces —49.13 MeV of the tensor force matrix 
element, which exhausts approximately 90% of the total value of —54.55 MeV in 
Table ILT1 This feature was first pointed out in ATMS. 25 ) In particular, three specific 
configurations, (0pi/ 2 )i , (lsi/2) (0^3/2)10 and (°P3/2) (0/5/2)10 for particle states, give 
large contributions in Table IVl They also have large mixing probabilities in the wave 
function in Table El These facts denote that three configurations are essential for 
describing the tensor correlation in 4 He. It is noted that the mixing probability of 
each configuration in Table IVTl is less than 3%, but their contributions to the tensor 
force matrix element (Vr) are relatively large. This feature was also discussed in the 
previous paper. 4 ) 

§5. Conclusions 

We have developed a method of calculating the nuclear ground state using the 
nucleon-nucleon interaction in the shell model framework. The important features of 
the nucleon-nucleon interaction are the strong tensor interaction caused by the pion 
exchange and the strong short-range repulsion caused by the internal structure of 
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Table VI. Contributions of each 0p0h-2p2h coupling to (Vr) in MeV. 



two particle states 


{$0p0h\V T \$2p2h) 


(Opi/ 2 )io 


-8.38 


(lsi/ 2 )(Od 3 / 2 )io 


-10.99 


(PP3/2) (O/5/2) 10 


-10.17 


(Orf 5 / 2 ) (Op 7 / 2 )io 


-5.67 


(0d 3 / 2 )? 


-2.62 


(Opi/ 2 )(Op 3/2 )io 


-2.48 



the nucleon. We have treated the tensor interaction in terms of the tensor-optimized 
shell model (TOSM) in which, in addition to the core state, we introduce two- 
particle two-hole (2p2h) states, to take into account the excitations due to the tensor 
interaction. As for the short-range repulsive interaction, we have introduced the 
unitary correlation operator method (UCOM), in which the unitary transformation 
is introduced to express the short-range behavior of the relative wave function. We 
have then combined these two methods to calculate the nuclear ground state, which 
is called TOSCOM. 

In TOSCOM, we have worked out the formulation in the Gaussian basis func- 
tion, and all the matrix elements are given in the Appendix. We obtained good 
convergence in the calculated results of the energy and other components of the 
Hamiltonian. This denotes that short-range correlation is successfully described us- 
ing UCOM, and simultaneously, the tensor correlation is also described explicitly 
using TOSM with the inclusion of the 2p2h states with high-momentum. Hence, 
we have developed the method of describing the nuclear structure starting from the 
nucleon-nucleon interaction in TOSCOM. 

We have carried out calculations first for 4 He using the UCOM for all the partial 
waves. We have found that the tensor correlations were largely underestimated and 
the binding energy of 4 He is also somewhat underestimated. The reason is the 
removal of the short-range part of the relative wave function in the optimization 
of the short-range correlations, where the tensor interaction needs some amount of 
strength. To overcome this feature, we have newly introduced S'-UCOM, where 
the unitary transformation was performed only for the s-wave component of the 
relative wave function. It was shown that the situation was largely improved, and 
the numerical results are found to be very close to the rigorous calculation. This is a 
very encouraging result to describe nuclei using the nucleon-nucleon interaction. We 
have, however, still somewhat lack of contribution of the tensor interaction matrix 
element. 

We consider that we have carried out the best calculation in the present frame- 
work. We have introduced one approximation of truncating the correlated operators 
due to the unitary transformation up to the two-body terms. This truncation has 
been shown to be good in the case of the central interaction alone by Feldmeier 
et al., 8 ) since the short-range repulsion is of short-range. However, the results of 
the tensor interaction obtained in the TOSCOM formulation seems to require more 
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attention because of the interference with the short-range repulsion. The tensor in- 
teraction seems to require some short-range components, which are taken away by 
the UCOM treatment of the short-range repulsion. We are currently investigating 
the effect of the three-body terms and the results will be reported in the near fu- 
ture. It could also be an idea to perform further UCOM with the tensor correlation 
operator 9 ^ in the very short-range part on top of TOSCOM, which can treat the 
intermediate- and long-range parts of the tensor correlations. We would like to note 
here that the present calculation is nearly a variational calculation in the shell model 
basis, and the numerical results are very encouraging for expressing all the neces- 
sary correlations in the calculated wave functions caused by the nucleon-nucleon 
interaction. 
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Appendix A 

Two-Body Matrix Elements in TOSM 

We write here the two-body matrix elements of the central, LS and tensor inter- 
actions in the Gaussian basis function. We also expand the potential with the finite 
number of the Gaussian function. Hence, we need to calculate the two-body matrix 
elements of the potential having Gaussian form. We do not write the isospin part. 
The matrix elements are calculated by transforming wave functions from jj coupling 
scheme to LS coupling scheme. In the following, we define L as the coupled orbital 
angular momentum of the two-particle states in the ^-coordinate bases. 

First, we define the Gaussian basis function for one nucleon state having the 
orbital angular momentum l,m and the length parameter b. 

ut m (r) = N l (a)r l e-% r2 Y lm (r), a = ±, N t (a) = 

Using this basis function, we evaluate the formulae of the matrix elements for the 
central, LS and tensor interactions and further the overlap between V and T-type 
bases for s-wave two-nucleon relative motion in the T-type bases. 

A.l. Central interaction 

We write the central interaction as 



2 a'+ 3 / 2 

r(i + 3/2) 



(Al) 
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We also consider multiplying the factor r^ry 2 , which is used for ls-wave state. The 
matrix element is given as 

4 



{ _ X) h + U-L H{ Nli(ai) .l^ 



i=l 



Y,W(hl2hh;Lk)(h0l 3 0\kO)(h0U0\kO) jf™ k (p,/3i,/? 2 ). (A-3) 



fc=0 



Here, W (hfohh; Lk) is Racah coefficient and ij^^'^^p, (3i, (3 2 ) is defined as 



Pi -A» 
2 



A»-Pi 
2 



V A 



2 \ n 



-^P2+A*+2+2n (A) (A-4) 



with 



Pi = Zi + Z 3 + ni, P 2 = / 2 + /4 + n 2 , / = \/2Z + 1, 

oi +q 3 q 2 + q 4 p 2 

Pi = z +P, /3 2 = z + P, P3 = 02--Z-: 

11 Pi 



where 



n 



n+l - 
2 - 



~ n+l 

2a~ 



^-i(a) 



2 ' 



Here, (x) n = x(x + l)(x + 2) • ■ ■ (x + n — 1) = n™=o + * 

and (x)o = 1. 

A. 2. LS interaction 

We write the LS interaction as 



r(x + n) 



(A-5) 
(A-6) 

(A-7) 
where x > 0, 



y L5 (r) = e -^i-^) 2 (x . 5). 



(A8) 



Here 



L = rxp, r = r 1 -r 2 , p = -(Pi~P 2 ), S = si + s 2 . (A-9) 
The matrix element of the LS interaction is given as 

e^ r2 {L-S) 

(-l)^i-J > /6W(LlL / l;Jl)<[< 1 ,«£'] L || e -^X||[«^,«J] I/ > J (A10) 
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where xi is the coupled wave function of two intrinsic spins with triplet state. The 
reduced matrix elements including the orbital angular momentum L consist of four 
terms as 



<K.««]L||e" 



L \\l u l 3 ' U l 4 \L> 



([O^lle-^ h\\K^u%] L ,) + ([ U %, U %} L \\e-^ Z 2 ||[^,<» 



-([ul\u%} L \\e-? r \ 2 xpMul 3 ,u%} L ,) 



(A-ll) 



Here, the first term including Zi is 



h ^ l2 ]L\\e- pr <iiik>?» 



r- 4 

= ^ ■ II • (-^"'^Vfeas + l)-^' • W(l 3 Ll 3 L';l 4 l) 



i=i 



(Mhfh ^(/i0/ 3 0|/x0)(Z 2 0/40|/x0)W(/i/2/3i4; L/i) 



^£L>/^ 2 ), (A-12) 
and the second term including Z 2 is 

(K 1 ,Mj 2 ]L||e~ pT ' 2 i2||[«j3,«k]L'> 

= ^ II ^ • (-l) Il_ ' 4+L+L '"V^(/4 + l)- L L> ■ W(hLl 4 L>; l 3 l) 
i=i 

x /ir 2 4(4) 2 ^(/i0/ 3 0|^0) (Z 2 0Z 4 0|/x0) WihhhU-Lp) 



X 



The other reduced matrix elements are 

(K 1 ,< 2 2 ] L ||e-^ 2 r 1 xp 2 
A" 4 

= ^ ■ II N k(°i) -VG-LL'-h h h £ (-i)^ A3+L <Z 3 010|A 3 0) 

1=1 A3,A4 

z 3 14 i' f • A 3 A 4 E(/i0A30|^0)(/ 2 0A40|/i0)M/(ZiZ 2 A 3 A4;L/i) 



(A-13) 



l 



° 4 ZfiMsjM fa, /h) (04 + 1 • <*A4,i4+i - ■ 8\t,u-i) 



+ y/h ■ (2l A + 1) l£l£julP> A) ■ 



(A-14) 
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and 



7T 



<K\< 2 2 ]L||e-^ r 2 x Pl ||[^,<» 

4 

H^K) -Ve-LL'-h i 2 f 4 (-i) h - X3+L (/ 4 oio|A 4 o) 

i=l A 3 ,A 4 
A3 A 4 L "| 

z 3 k z/ }> ■ A 3 A 4 ^(/i0A30|/x0)(Z20A40|/i0)W(ZiZ2A 3 A4;-LAi) 



1 1 1 



°3 ^Jkk^ 1 '^ 2 ) (03 + 1 • <Ja 3 ,/ 3 +1 - V^3 • «5a 3 ,/ 3 -i) 

+ • (2z 3 + 1) i?;j 2 )lM h, A) • Wi 



Here, {• • • } including nine numbers is 9j symbol. 

A.3. Tensor interaction 

We write the tensor interaction as 

where 

3(<ri • r)(er 2 • r) 



S\2 — 



— <J\ ■ cr 2 



24vr 



The matrix element of the tensor interaction is given as 



K,u%] L , X i 



:< 3 3 ,<>,xi 



where 



r m e - P v Su 

= 4 V^{-l) L+1 ~ J W(L1L'1; J2) 

x {[u h ,u l2 ] L \\r m e- pr2 Y2(r)\\[u h ,u u ] LI ), 



(K,v%] L \\r m e-" y 2 (f)||[< 3 3 ,<>) 
4 



P1,P2, X 1< L 1 
«l='l-Pl. 
12= l 2-P2 



P 3 >P4. A 2> i 2 
9 3 =i 3 "P3. 



x (_i)A 2 +l 2 X1A2L1L2 ■ J]<AiOA 2 0|aO) s 

S fj, 

( Ai Li L ^ 
x (Li0L 2 0|^0)(^0s0|20)<^ A 2 L 2 L' i 



4-r(/i + |) 



Pl+P 2 ,<?l+9 2 ,P3+P4,<?3+<?4 I 2' a '° 2 J ' 



S /i 2 





(A-15) 



(A-16) 



([cri,cr 2 ] 2 .y 2 (f)). (A-17) 



(A-18) 



(A-19) 
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Here, p 1 + q± = l\ and p 2 + <?2 = k, and 

c i uh ,L _ (~1) P2 / (2Z 1 + 1)!(2Z 2 + 1)! 

PUP2MM 2P1+P2 y ( 2pi + i)!(2 9l + l)!(2p2 + 1)K 2 92 + 1)! 

„ „ f pi 9i ^1 
x pufegig 2 iii2< P2 92 la ^ (pi0p 2 0|Ai0)(gi0g 2 0|Li0), (A-20) 
L Ai Li L J 



a = + a 2 + a 3 + a 4 ), /3 = ^(-ai + a 2 - a 3 + a 4 ), (A-21) 

a 1 = - +P , a2 = a -— I -. (A-22) 

A.4. Overlap between V -type basis and T-type basis 

We consider the overlap between V-type basis and T-type basis. Here, we limit 
the case with s-wave relative motion for the T-type basis, whose length parameters 
are a r and or for r, R of the T-type basis, respectively. The overlap matrix element 
is given as 

L oo , — 

= N h ( ai ) N l2 (a 2 ) N (a r ) N L (a R ) £ ^ (- 1 )' 2 " 9 

g=0 A=0 

x W(h k qL-q;LX) l\kq(L- q) (Z^AO) l\^- q (f . «> , (A-23) 



where 



a = flr + T' ^ = ar_ T' A = Zi+g + ni, (A-24) 



,l 1 (2L + 1)! 



P 2 = Z 2 + L-, + n 2 , = Y^ (2q + l)\.(2L-2q + l V : < A " 25 > 

The overlap with ci-wave relative motion for r in T-type basis can be calculated in 
the same manner. 
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